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ABSTRACT 



Aims. The recent discovery of the radio shell-type supernova remnant (SNR), G353. 6-0.7, in spatial coincidence with the unidentified TeV source 
HESS J173 1-347 has motivated further observations of the source with the High Energy Stereoscopic System (H.E.S.S.) Cherenkov telescope 
array to test a possible association of the y-ray emission with the SNR. 
^ Methods. With a total of 59 hours of observation, representing about four times the initial exposure available in the discovery paper of 
^ HESS J173 1-347, the y-ray morphology is investigated and compared with the radio morphology. An estimate of the distance is derived by 
^ comparing the interstellar absorption derived from X-rays and the one obtained from ^^CO and HI observations. 

^ Results. The deeper y-ray observation of the source has revealed a large shell-type structure with similar position and extension (r~0.25°) as the 
radio SNR, thus confirming their association. By accounting for the H.E.S.S. angular resolution and projection eff'ects within a simple shell model, 
the radial profile is compatible with a thin, spatially unresolved, rim. Together with RX J1713.7-3946, RX J0852.0-4622 and SN 1006, HESS 
J173 1-347 is now the fourth SNR with a significant shell morphology at TeV energies. The derived lower limit on the distance of the SNR of 
3.2 kpc is used together with radio and X-ray data to discuss the possible origin of the y-ray emission, either via inverse Compton scattering of 
electrons or the decay of neutral pions resulting from proton-proton interaction. 
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1. Introduction 

In the survey of the Galactic plane carried out by the H.E.S.S. ex- 
periment, many sou rces emitting at T eV energies remain uniden- 
tified to date (e.g. lAharonian et alJ 12008). Most of the sources 
are extended beyond the point spread function (PSF) of the 
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H.E.S.S. experiment (~ 0.06° for the analysis presented in this 
paper). The largest number of conclusive identifications so far 
can be attributed to pulsar wind nebulae (PWNe) as presented 
in e.g. iGallant et aP fcOOSh . Recently, a ne w radio SNR, ca ta- 
logued as G353.6-0.7, was discovered by .Tian et al.l (1200 8^ to 
be in spatial coincidence with H ESS J173 1-347 , one o f the 
unidentified sources presented in Aharonian et al.l (l2008h . The 
diameter of the radio shell is nearly 0.5° which allows, given 
the brightness of the source and the H.E.S.S. angular resolu- 
tion of -0.06°, for a morphological comparison of the y-ray 
source with the shell observed in radio. Moreover, at least up 
to the current date, no radio pulsar or X-ray PWN candidate 
was found that might alternatively explain the TeV emission. 
This situation should be compared to other y-ray sources lik e 
HESS J1813-178 or HESS J1640-465 (lAharonian et al.ll200d) . 
which are also in spatial coincidence with radio SNR shells. 
However, for these latter sources a morphological identifica- 
tion with the radio shells is not possible, and the emission can 
pl ausibly orig i nate fr om a PWN seen in X-rays as discussed 



in 



Funk et ai ; ('2007b') for HESS J1640-465 and in i Funk et al.1 



(l2007ah :[Gotthelf & Halpern (2009) for HESS J1813-178. 

Observations of the north-eastern part of HESS J173 1-347 
with the X-ray satellites XMM-Newton, Chandra, and Suzaku 
have confirmed an X-ray counterpart found in archival RO SAT 
data (presented in 'Aharonian et al.ll2008l: iTian et al.ll20Q8 ). An 
X-ray shell partly matching the radio morphology was found and 
the spectral analysis has revealed that the X-ray emission is of 
synchrotron origin, indicating that the shoc k wave of the SNR 
has accelerated e lectrons up to TeV energies (lAcero et al.l2009bl: 
iTian et all l2010). A compact (unresolved) X-ray source XMMU 
J173203.3-344518 (Halpern & Gotthelf 2010) was observed to- 
wards the geometrical center of the remnant and has spec- 
tral properties reminiscent of central comp act objects (C CQs) 
found in several other supernova shells (e.g. Pavlov et"aD l2004). 
A search for pulsations using the EPIC PN cameras onboard 
XMM-Newton shows on ly marginal evidence of a 1 s period 
(iHalpem & Gotthelfll2QrO ). 

Given the recent disco very of G3 5 3. 6-0. 7, little is known 
about its age and distance. iTian et all (l2008l) suggested a dis- 
tance of 3.2 ± 0.8 kpc assuming that the SNR is at the same 
distance as the Hii region G353.42-0.37. 

Additional H.E.S.S. observation s, carried out since the dis- 
covery paper of HESS J1731-347 (lAharonian et al.l 12008 '). al- 
low to investigate the compatibility of the TeV source with the 
radio shell SNR G353.6-0.7. The observations and the data anal- 
ysis are described in Sect. |2] and the morphological and spectral 
results in Sect. [3] The multi-wavelength counterparts of HESS 
J173 1-347 are described in Sect. |4] and a general discussion is 
presented in Sect. [3 



2. H.E.S.S. observations and analysis methods 

H.E.S.S. is an array of four identical imaging atmospheric 
Cherenkov telescopes (lACTs) locat ed in the Khomas Hig hland 
of Namibia 1800 m above sea level (iBemlohr et al.ll2003l) . The 
survey of the Galactic plane by the H.E.S.S. collaboration has 
led to the discovery of the y-ray source HESS J173 1-347, pre- 
sented as an unidentified extended source in lAharonian et al.l 
(l2008l) . In this first data set, 14 hours of observation time were 
available. Additional dedicated observations were carried out 
in July 2007 and in July and August 2009 with zenith angles 
ranging from 9° to 42°, the mean angle being 16.5°. The total 
H.E.S.S. observation time for this target is 59 hours after data 
quality cuts. 



The data set was ana lyzed using the Model analysis 
(Ide Naurois & Rollandll2009 ^ which exploits the full pixel infor- 
mation by comparing the recorded shower images with a pre- 
calculated shower model using log-likelihood minimization. In 
comparison with conventional analysis techniques, no cleaning 
or parametrization of the image shape is required and the full 
camera information is used. This method leads to a more accu- 
rate reconstruction and better background suppression than more 
conventional techniques and thus to an improved sensitivity. 

Spectral and spatial analyses were carried out using a 
minimum image intensity of 60 photoelectrons (p.e.) resulting 
in an energy threshold of 240 GeV and an angular resolution 
of 0.06° (68% containment radius). All resul ts presented were 
cross-checked with a multivariate analysis (lOhm et al.l l2009l) 
using an independent calibration and gamma/hadron separation, 
which yielded consistent results. Unless otherwise quoted, the 
error bars in the following section are given at 1 cr. 



3. TeV y-rays analysis results 

The H.E.S.S. excess map of the region of HESS J173 1-347 
is shown in Fig. [T] smoothed with a Gaussian of (T=0.04°. 
For the background estimation in the image and in the mor- 
phology studies, the ring background method presented in 
Berge et al.' ("2007") was used. Because of the larger data set and 
the more sensitive reconstruction technique, the presented im- 
age is mu ch more detailed than the one shown in the discov- 
ery paper (lAharonian et al ]l2008l) . This reveals a complex region 
composed of a large and bright structure (HESS J173 1-347), 
detected at 22(T, with a suggestive shell-like morphology. 

A smaller and fainter structure named HESS J1729-345 
(detected at 8(t) is also observed, the properties of which are 
presented separately in Sect. 13.31 



3.1. TeV energy morphology 

To further test the hypothesis of a shell morphology for HESS 
Jl 73 1-347 and its association with the radio SNR, radial and 
azimuthal profiles in radio and y-rays were extracted centered 
on the position of the CCO (aj2ooo =17h32m03s, ^j2ooo = 
-34°45'18'0, also coincident with the geometrical center of the 
radio SNR. The profiles were derived from the uncorrelated y- 
ray excess map and corrected for the field of view (FoV) accep- 
tance. 

For the y-ray radial profile, the position angle£] from 270° 
to 310° were excluded, to avoid contamination from HESS 
J1729-345, and the resulting radial profile was compared with a 
sphere and a shell model. The first model is a uniformly emitting 
sphere of adjustable radius, projected on the sky and then folded 
with the PSF derived for this analysis (r68% = 0.06°). The shell 
model consists of a uniformly emitting shell of variable outer ra- 
dius and thickness (defined as router - Hnner) projected on the sky 
and then folded with the same PSF. 

In the morphological test, the best fit statistically favors the 
shell model and the sphere model is ruled out at 3.9 cr {x^l&of = 
2.90/5 and 28.12/6 for the shell and sphere models respectively). 
In the case of a shell model, the best fit radius is 0.27° ± 0.02° 
and the emission is compatible with a thin, spatially unresolved, 
shell with an upper limit thickness of 0.12° (90% confidence 
level). 



Position angle 0° corresponds to North and 90° to East. 
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Fig.l. TeV r-ray excess map (1.5° x 1.5°) of the HESS 
J173 1-347 region smoothed with a Gaussian width (T=0.04°. 
The average H.E.S.S. PSF for the dataset is shown in the inset. 
The regions used for the spectral analysis of HESS J173 1-347 
and HESS J1729-345 are respectively represented by the large 
and small dashed circles. The position of the central compact 
object detected in X-rays is shown with a white cross. The linear 
scale is in units of excess counts per smoothing Gaussian width. 
The transition between blue and red in the color scale is at the 
level of 4cr. 



To compare the TeV morphology with the shell seen in ra- 
dio, the radio continuu m map from the ATCA southern Galactic 
plane survey (SGPS) (iHaverkom et al. 1 12006 ) was smoothed to 
match the H.E.S.S. spatial resolution and a radial profile was 
extracted (excluding point sources). The radio profile was then 
scaled by a normalization factor calculated as the ratio of the 
total number of excess y-rays over the total radio flux on the 
whole remnant. The resulting profiles, presented in Fig. |2] show 
an extended emission in y-rays similar to that seen in radio. 

In contrast with RX J17 13.7-3946 which is brighter in the 
North- West and SN 1006 that exhibits a bipolar morphology, the 
azimuthal profile of HESS J173 1-347 (see Fig. [3]) integrated for 
r < 0.3° shows no significant deviation from a flat profile (x^/dof 
= 8.8/9). 



3.2. Spectral results 

The energy spectrum of the SNR was o btained by means of a 
forward-folding maximum likelihood fit (iPiron et al.ll2001) from 
a circular region of 0.3° centered on the CCO, illustrated by the 
large dashed circle (r = 0.3°) in Fig.[Tl chosen to fully enclose 
the emission of the remnant. The background is estimated us- 
ing the multiple reflected-regions technique where background 
events are selected from regions of the same size and shape as 
the source reg ion and at equal an gular distance from the observa- 
tion position teerge et al.ll2007h . The resulting spectrum, shown 
in Fig. m is well described by a power-law model (equivalent 
X^ldof = 27.7 / 35) defined as dN/dE = No(E/Eo)-^ where Eq 
is the decorrelation energy (energy at which the correlation be- 
tween the slope and the normalization vanishes). The best fit pa- 
rameters, listed in Table [T] result in an integrated 1-10 TeV en- 



ergy flux of (6.91 ± 0.75stat ± 1.38 



syst 



) X 10 erg cm ^s 
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fl ux measured here is low er than what has been derived initially 
in lAharonian et al ] (l2008h : (16.2 ± 3.6stat ± 3.2sys,) x IQ-'^ erg 
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Fig. 2. The y-ray excess and radio radial profiles are shown 
with green crosses and red squares respectively. The best fits 
to the y-ray data of a sphere and a shell model are overlaid. 
Both radial profiles are centered on the compact central object 
(c^j20oo =17h32m03s, ^j2ooo = -34°45q8'0. 
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Fig. 3. Normalized azimuthal y-ray excess profile restricted to 
radius r < 0.3° and using the same center as in Fig. [21 The bright- 
ness distribution is compatible with a flat profile. 



cm"^s"^ in the same energy band. However, the region of ex- 
traction in the discovery paper was much larger (r = 0.6° versus 
r = 0.3° in this paper), including HESS J1729-345 and pos- 
sibly some surrounding diff'use emission. A cross-check to de- 
riv e the flux from th e SNR only using the same data set as used 
in lAharonian et al.l (l2008k and following the original analysis 
method gave results consistent with the complete data set pre- 
sented here thus confirming that the flux difl'erence was mainly 
due to the choice of the integration region. A power-law model 
with an exponential cutofl' was also tested which did not improve 
the quality of the fit (equivalent ;^Vdof = 24.0 / 34). 



3.3. HESS J1 729-345 

A y-ray excess of TeV emission was found at the best fit position 
o^moo =17h29m35s, feooo = -34°32'22'' with a statistical error 
of 0.035°and the source was therefore labeled HESS J1729-345. 

The source is extended beyond the size of the PSF (Gaussian 
width cr = 0.12° ± 0.03°) and the region used to derive the spec- 
tral parameters is shown by the small dashed circle (r = 0.14°) 
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Fig. 5. Multi- wavelength view of the HESS J173 1-347 region. The radio and the y-ray image show the same field of view while the 
X-ray image is zoomed in order to show the details of the shell structure. The significance contours at 4, 6 and 8 cr obtained with an 
integration radius of 0.06° and the Galactic plane (white dashed line) are overlaid in the three panels. Left : ATCA radio map at 1.4 
GHz from the south Galactic plane survey (SGPS) in units of Jy/beam with a beam of 100". The HII regions G353.42-0.37 (left) 
and G353.38 1-0.1 14 (right) are marked with arrows. Middle : XMM-Newton observation of a sub-region of the SNR, in the 0.5-4.5 
keV energy band, using MOS instruments with units in ph/cm^/s/arcmin^. The position of the source XMMU J173203.3-344518, 
which is likely to be the CCO of the SNR is shown by the red arrow. Right : TeV y-ray excess map of HESS J173 1-347 smoothed 
with a Gaussian with cr=0.04°. The region used to derive the radio flux and the spectral parameters in X- and y-rays for the SED is 
also shown. 
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Fig. 4. Difl'erential energy spectra of HESS J173 1-347 (filled 
circles) and HESS J 1729 -345 (open circles). The normalization 
for the second source has been divided by 10 for graphical pur- 
poses. Events were binned to reach a significance of at least 2cr. 
The best fit power-law models along with the residuals for HESS 
J173 1-347 are also shown. The grey bands correspond to the 
range of the power-law fit, taking into account statistical errors. 

in Fig. [U The spectrum obtained is well modeled by a power- 
law model (see Fig. (4]) and the best fit parameters are listed in 
Table [T] The integrated flux in the 1-10 TeV energy band is 
(0.88 ± 0.29stat ±0.18syst) X lO'^^ erg cm'^s'^. 

4. Multi-wavelength counterparts 

One of the interesting characteristics of HESS J1731 -347 is that 
non-thermal emission is clearly identified in radio. X-rays and 
at TeV energies. In X-rays however, the access to the spectral 
properties is limited to a subregion of the SNR as the coverage 
with the XMM-Newton, Chandra and Suzaku satellites is only 



partial, and the statistics in the ROSAT All Sky Survey data are 
too low. In order to study the spectral energy distribution (SED) 
of the source, the radio flux and the TeV spectral properties were 
extracted only from the region observed in X-rays (see region 
definition in Fig. [51 right). The multi- wavelength counterparts of 
HESS J1729-345 are discussed later in Sect.O 

4.1. Radio Continuum 

The shell observed in radio is spatially coincident with the y-ray 
shell and has a similar extent (see radial profile in Fig. O. The 
flux obtained (excluding point sources) from the SGPS ATCA 
data in the region observed by the XMM-Newton pointing is 0.8 
± 0.3 Jy at 1420 MHz. The total radio flux for the SNR measured 
bv lTianetal.1 (12008 ) is of 2.2 ± 0.9 Jy. The compact HII region 
(G353.42-0.37) l ocated to the We st of the remnant at a distance 
of 3.2 ± 0.8 kpc (iTian et al.ll2008h is indicated in Fig. [5] (left). 

4.2. X-rays 

In order to derive spectral information from the X-ray emis- 
sion from the remnant, the XMM-Newton pointing obtained as 
a follow up of the HESS source (Obsid: 0405680201 ; PI: G. 
Ptihlhofer) was analyzed. To clean the proton flare contamina- 
tion during the observation, a histogram of the 10-12 keV count 
rates of each camera was built. A Gaussian fit was then per- 
formed in order to remove time interval s where the count rate s 
were beyond 3 cr from the mean value (iPratt & Amaudll2002h . 
The remaining exposure time after flare screening is 22 ks out of 
the 25 ks of observation for MOS and 15 ks for PN. For the im- 
age generation, the instrumental background was deri ved from 
the compilation of blank sky observations by Carter & ReadI 
( 2007) and renormalized in the 10-12 keV band over the whole 
FoV The image resulting from the combination of the two MOS 
instruments is presented in Fig. [5] (middle). For this mosaic, the 
data from the PN instrument were not used because of stray light 
contamination to the North-East (photons singly reflected by the 
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Table 1. Best fit spectral parameters obtained for diff'erent extraction regions in HESS J173 1-347. The model used is a power-law 
of the form dN/dE = No(E / Eo)~^ . The systematic errors are conservatively estimated to be ± 0.2 on the photon index and 20% on 
the flux. 



Region 


Photon index F 


Decorrelation energy Eq 


Normalization A^o 


1-10 TeV integrated flux 






TeV 


10-12 cm-^s-i TeV-i 


10-12 erg cm-^s-i 


HESS J1731-347 


2.32 + 0.06stat 


0.783 


4.67±0.19stat 


6.91 ± 0.75stat 


sub-region of HESS 11731-347^ 


2.34 + 0.09stat 


0.780 


1.41±0.11stat 


2.02 + 0.36stat 


HESS J1729-345 


2.24±0.15stat 


0.861 


0.44 ± 0.07stat 


0.88 + 0.29stat 



^ A spectral analysis corresponding to the FoV of the XMM-Newton data (see Fig.[5l center) has been carried out in order to build a SED. 



mirrors) from a bright X-ray source located outside the FoV. This 
results in some spurious arc features near the border of the FoV 
in the North-East. 

The X-ray emission is characterized by extended emission 
which is concentrated in arc-like features, similar to broken 
shell seen from many shell-type SNRs. Some of the arcs partly 
coincide with the radio and y-ray shell (see Fig. O. Some of 
the structures could hint at an additional, smaller shell, but 
might also come from irregular SNR expansion in an inhomo- 
geneous and/or dense medium (iBlondin et al.ll200i1) . A double- 
shell structure is also observe d in RX Jl 7 13. 7 -3946 in X-rays 
(iLaz endic et al. 2004; Cassam -Chenai et al ] 120041: lAceroet al.l 
[2009ah . 

The spectral analysis of the difl'use X-ray emission was car- 
ried out using the Extended Source Analysis Software (ESASQ) 
provided in the XMM-Newton Science Analysis System (SAS 
v9.0) to model the particle and instrumental backgrounds. The 
error bars in this section are quoted at 90% level confidence. 
For this analysis, the three instruments PN+M0S1+M0S2 were 
used and the regions were selected to avoid the stray light fea- 
tures. 

The spectrum derived from the region covered by the FoV of 
XMM-Newton that is used for the SED is shown in Fig.[6j The 
emission is well represented by an absor bed power-law m odel 
and no emission lines were found (see alsol Tian et al.ll2010l) . The 
best fit parameters obtained from a joint fit of MOSl, M0S2 and 
PN spectra are A^h = (1.08 ± 0.02) x 10^^ cm"^, a spectral index 
F = 2.28+0.03 and a normalization at 1 keVA^o = (1.37±0.05)x 
10"^ cm"^ s"^ keV"^ A search for spatial spectral variations of 
the diff'use emission revealed that the power-law index is in most 
locations in the range F = 2.1 - 2.5. Under the assumption of a 
pure power-law hypothesis, the absorption column significantly 
increases towards the Galactic plane from A^h = (0.93 ± 0.05) x 
10^2 cm-2 in the South-East region to A^h = (2.23 ± 0.21) x 10^^ 
cm"^ in the North- West region (see Fig. [7]; left). The errors on 
the absorption column shown in Fig. [7] (left) are ranging from 
5% to 12%. 

The bright point source XMMU J173203.3-344518 lies at 
the geometrical center of the radio and y-ray shell. Marginal ev- 
idence for a pulsation at a period of 1 s for the pulsar candi- 
date and a faint nebula (radius of 30") whose spectral properties 
are compatible with a dust - scattered halo have been reported by 
^Halpem & Gotthell (12010 ). As no optical or IR counterpart of 
the point source have been detected and as the X-ray spectrum 
is well described by a blackbody emission model with kT ~ 0.5 
keV, the obj ect is a g ood candidate to be the CCO of the SNR 
(lAcero et al . 2009b; Halpem & Gotthelfll2010l:lTian et al.ll2QlQl) . 



^ http : //xmm2 . esac . esa. int/external/xmm_ s w_cal/background/ 
epic_esas.shtml 
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Fig. 6. X-ray spectrum, using the PN camera, extracted from 
the sub-region of HESS J173 1-347 shown in Fig. [5] (right). 
The non-thermal emission from the SNR is described by an ab- 
sorbed power-law model (dashed line). The local astrophysical 
background, fitted to an off' region outside the SNR, is mod- 
eled by two components (dotted lines). The low energy com- 
ponent is an APEC model (astrophysical plasma emission code, 
see http://hea-www.harvard.edu/APEC) representing the back- 
ground from the Local Bubble and the high-energy component is 
an absorbed power-law representing the hard X-ray background 
(unresolved AGNs, cataclysmic variables, etc). The residuals 
of the total model (SNR+local astrophysical background) are 
shown in the lower panel and the;^Vndof is 1921/1569. 



4.3. ^^CO(J=1-0)andHI 

The comparison of the absorption along the line of sight de- 
rived from X-ray data, ^^CO and HI observations can be used 
to constrain the distance to the SNR. The velocity s pectra of the 
^^CO emission (using data from the CfA survey, iDame et aP 
200 ID and the HI emi ssion (using data from the SGPS survey, 
Haverkom et al. derived from the region of highest X- 

ray absorption (aj2ooo =17h31m43s, ^j2ooo = -34°34'58'0 are 
shown in Fig. [8] (bottom). 

In order to derive a lower limit on the integration distance 
required to match the A^h derived from X-rays, all the material 
is assumed to be at the near distance allowed by the Galactic 
rotation curve. Under this hypothesis, the cumulative absorption 
column derived from the atomic and molecular hydrogen shown 
in Fig. [8] (top) is similar to the one observed in X-rays, Nu = 
(2.23 ±0.21) X 10^^ cm"^, when integrating up to a radial velocity 
relative to the local standard of rest (LSR) of -25 km/s. The CO- 
to-H2 mass conversion factor and the HI brightness temperature 
to column density used are respectively of 1.8 x 10^^ cm"^ K"^ 
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Fig. 7. Left : X-ray absorption map derived from a spectral fit to XMM-Newton data assuming an absorbed power-law model. 
A significant increase of A^h towards the Galactic plane is observed. Right : Absorption column map derived from atomic and 
molecular hydrogen when integrating over radial velocities from km/s to -25 km/s (see Sect. l4.3l for more details). The Galactic 
plane is represented by the white dashed line. In both panels, the XMM-Newton field of view is represented by a dashed circle and 
the X-ray contours obtained from Fig. (center) are overlaid. 



km"^ s ("Dame et al."200l') and 1.82 x 10^^ cm'^ K'^ km'^ s 
(iDicke v & Lockman 1990). 

When integrating up to the same velocity, the map of Nw 
derived from the atomic and molecular hydrogen shown in Fig. [7] 
(right) exhibits an increase of absorption towards the Galactic 
plane similar to that in the X-ray absorption map in Fig. |7](left). 
The peak of ^^CO emission at a LSR velocity of -18 km/s is 
thus in the foreground of the SNR and is likely to be the cause 
of the gradient of absorpti on seen in X -rays. Using the circular 
Galactic rotation model of lHou etal.1 (12009) with a distance to 
the Galactic center of 8.0 kpc, the nearest distance corresponding 
to the LSR velocity of -18 km/s is 3.2 kpc thus setting a lower 
limit for the distance of the remnant. 

4.4. GeVy-rays 

In the Fermi-LAT first year catalog (lAbdo et al.ll2010h the source 
IFGL J1729.1 -3452c is found in the neighborhood of HESS 
J173 1-347 as shown in Fig.[9l The Fermi source has an analysis 
flag that indicates that the source position moved beyond its 95% 
error ellipse when changing the model of difl'use emission. The 
^^CO map in Fig. [9] shows that the Fermi source is located near 
a small scale gas clump that could be not well represented in the 
difl'use emission model. The position of the source presented in 
the catalog is to be used with caution and is therefore possibly 
not incompatible with HESS J1729-345. 

The Fermi source has a photon spectral slope of 2.26+0.08 
and shows neither indication for spectral curvature nor time vari- 
ability on a time scale of months (the catalog does not address 
shorter or longer time variations). This source is the closest 
Fermi detection near the newly discovered SNR and the flux de- 
rived in the Fermi catalog is used as an upper limit in the SED 
of the SNR in Fig.IH 

4.5. Multi-wavelength counterparts for HESS J1 729-345 

At radio wavelengths, the y-ray contours of HESS J1729-345 
lie near the HII region G353.381-0.114. Using HI radio recom- 



bination line data, the LSR velocity corresponding to this source 
is either -54 km/s or -82 km/s (Caswell & Haynes 1987). In the 
latter case this HII region could be associated with the molecular 

cloud observed around velocities of 80 km/s (see Fig. [9]). At 

X-ray energies, no archival dedicated observations were found, 
and no emission is detected in the ROSAT all sky survey, proba- 
bly due to the high absorption in the line of sight. As discussed 
in the previous section, a Fermi source is found to lie close to 
HESS J1729-345. 



5. Discussion 

The newly discovered SNR HESS J173 1-347 is in several ways 
comparable to RX J1713.7-3946 and RX J0852.0-4622. Those 
objects are X-ray synchrotron emitters and exhibit no thermal 
emission lines. A CCO is also found within those three SNRs 
indicating a core collapse SN. Moreover at a distance of 3.2 kpc 
(see Sect.O, the TeV luminosity of HESS J1731-347 in the 
1-30 TeV energy band is 1.07x((i/3.2kpc)2 10^^ ^-i ^^j^h 
is similar to the luminosity of RX J17 13.7-3946 (the brightest 
TeV shell SNR detecte d until now), of 0.81x10^^ erg s"^ us- 
ing a distance of 1 kpc (Fukui et al. 2003:'Cassam-Che nai et al.l 
2004j; iMoriguchi et al.1 12005) and slightly higher than that of 
RX J0852.0-4622 with .65x10^^ erg s"^ at a distance of 0.75 
kpc (iKatsuda et al.ll200 8'). A difl'erence with RX J1713.7-3946 
is that the flat y-ray azimuthal profile of HESS J173 1-347 (see 
Fig. [3]) suggests that the remnant is evolving in a relatively uni- 
form ambient medium and that it is not in interaction with the 
cloud (shown in Fig. [9l top) used to derive a lower limit to the 
distance of the SNR. This significantly diff'ers from the case of 
RX J1713. 7-3946 which exhibits much brighter y-ray emission 
in the North- West where the shock is thought to interact with 
denser material. 

The distance used for the luminosity is derived from the ab- 
sorption in the foreground and provides only a lower limit of 
3.2 kpc. However, as it is believed that supernova explosions are 
more likely to occur in the spiral arms of the Galaxy w here the 
density of massive stars (i.e. SNR progenitors) is higher (iRusseill 
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Velocity (LSR) km/s 

Fig. 8. Top : Cumulative absorbing column density (solid line) 
as a function of radial velocity at the position of highest X-ray 
absorption (see Sect. 14.31) . The relative contributions from the 
atomic and molecular hydrogen are represented by the dashed 
and dash-dotted lines respectively. Middle : Rotation curve to- 
wards the same direct i on as derived from the model of Galactic 
rotation of Hou et al.1 (l2QQ9l) . Bottom : ^^CO (dashed line) and 
HI (dash-dotted line) spectra obtained the region highest X-ray 
absorption. 



l2003l:lHou et"aDl2Q09 '). it is likely that HESS J1731-347 could 
be located within the Scutum-Crux or Norma arms, which cross 
the lin e of sight at / = 353.5° at ~3.0 and ~4.5 kpc respec- 
tively (IHou et al.ll2QQ9l) . The next arm in the same line of sight 
is the Sagittarius arm lying at a distance of 12 kpc. This lat- 
ter possibility for the location of the SNR would lead to a 
much higher y-ray luminosity, an order of magnitude higher than 
RX J17 13.7-3946. Also at such a distance, the physical size of 
the remnant would exceed 50 pc, substantially larger than other 
TeV shell SNRs whose physical size is < 15 pc. As a result it 
is reasonable to believe that the real distance to the SNR should 
not be much larger than the derived lower limit of 3.2 kpc. 

The radio flux and the X- and y-ray spectra derived in Sect. 
131 from the sub-region of HESS J173 1-347 that is covered by 
the EoV of XMM-Newton were combined in the SED presented 
in Fig.[TOl The X-ray data were corrected for the interstellar ab- 
sorption with A^H = 1.08 X 10^^ cm"^. To mo del the SEP, a sim- 
ple one-zone stationary model (presented in lAcero et aL l l2QTQh 
was used. In this model, the spectrum of electrons and protons is 
represented by a power-law of slope s with exponential cutofl's 
at energies E^^q and £'c,p for the electrons and protons respec- 
tively. For the modeling of the object, it is assumed that the mea- 




17h32m 17h30m 17h28m 
Right ascension 




17li32m 17li30m 17li28m 
Right ascension 



Fig. 9. i^CO map of the vicinity of HESS J173 1-347 integrated 
from LSR velocity -13 km/s to -25 km/s (top) and from -75 
km/s to -87 km/s (bottom) respectively corresponding to the 
intervals of the second and third ^^CO peak shown in Fig. [S] 
The HESS significance contours from Fig. [5] together with the 
Galactic plane and the Fermi 95% position confidence level con- 
tours presented in the 1 year catalog are overlaid. The linear 
scale is in units of K km/s. 



sured multi- wavelength emission from the sub-region of HESS 
J173 1-347 is entirely coming from the SNR located at a dis- 
tance of 3.2 kpc. As this distance is only a lower limit, the total 
energy of accelerated particles (We and Wp) in the SNR should 
also be viewed as lower limits. 

In the pure leptonic scenario, the slope of the electrons is 
constrained by the radio and the X-ray synchrotron emission be- 
tween 1.9 and 2.1 and the strength of the magnetic field required 
to reproduce the ratio of observed synchrotron and IC emission 
Hes between 20 and 30 jdG for 15 < ^c,e < 25 TeV. Although 
the relative ratio of radio, X- and TeV fluxes can be fairly well 
reproduced by this leptonic scenario, the model is inadequate to 
account for the X-ray and the y-ray spectral slope as illustrated in 
Fig. [To] (top). The corresponding parameters for the latter model 
are summarized in Table |2] 
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This limitation no longer occurs in a scenario where the TeV 
emission is dominated by hadronic processes as the X- and y-ray 
emission are now independent and both spectral slopes can be re- 
produced as shown in Fig. [TOl (bottom). Moreover, the strength 
of the magnetic field can be increased as it is no longer fixed by 
the X/y ratio. In order to reproduce the observed TeV flux, the 
total energy in high-energy protons (E > I GeV) assuming a 
spectral slope of 2.0 is Wp = 2 x 10^^ (n/l cm-^)-^ (d/32 kpc)^ 
erg. It should be noted that this energy content only represents a 
sub-region of the SNR accounting for ~ 1/3 of the total TeV flux 
(see Table [TJ implying that the total energy transferred to accel- 
erated protons in the whole SNR is a substantial fraction of the 
energy available in the remnant for n ^ 1 cm"^ . For this energetic 
reason, gas densities much below this value appear incompatible 
with the hadronic emission scenario. 

Although it is not possible to measure the density of the 
ambient medium surrounding the SNR as no X-ray thermal 
emission is detected, an upper limit on the density can be de- 
rived. In order to do so, a thermal component, whose normal- 
ization is fixed for a given density using the method presented 
in Acero etal ] (12007) (Sect. 3.1), is added to the X-ray spec- 
trum. The shocked ambient medium is assumed to be in a non 
equilibrium ionization state with an ionization timescale param- 
eter T = 10^ cm"^ s and an electron plasma temperature kT^^l 
keV. Such values are commonly observed in other young SNRs 
for which the X-ray emission of the sho cked ambie nt medium 
has been studied as in e.g. RCW86 fsee lVink et al.l i2006). For 
the given parameters, the derived upper limit (90% confidence 
level) on the ambient medium density is 10"^ cm"^. In the case 
of a lower temperature (^Je^O. 15 keV), an upper limit of 1 cm"^ 
is reached. 

For a density of 1 cm"^ the corresponding shock speed and 
age of the SNR would be ~410 km/s and 14000 yrs in order 
to match a physical radius of 7?shock=15 pc (0.27° at 3.2 kpc), 
for a SN explosion of £^sn = 1 x 10^^ erg with a mass o f 
ejecta of 5 M© using equations from lTruelove & McKed (Il999l) . 
However, this shock speed is an order of magnitude lower than 
what has been measured in other bright synchrotron emitting 
SNRs like SN 1006 (Vsh^SOOO ± 400 km/s at a distance of 2.2 
kpc ; iK atsuda et al. 2009), RCW 86 (ysh=6 Q00 + 3000 km/s ; 
Helder et al...20091 CasA {V,^=A9QQ km/s ; iPatnaude & FesenI 
20091) or Tycho {V,^ = 30 00 ± 1000 km/s at a distance of 2.3 
kpc ; iKatsuda et al.ll201Qh . As a rough estimate, the required 
density to reproduce a canonical shock speed of 3000 km/s us- 
ing the aforementioned SN parameters is of the order of 0.01 
cm"^ (compatible with the upper limit derived from the lack of 
thermal X-ray emission in the previous paragraph) for a corre- 
sponding age of ~ 2500 yrs. 

To summarize, the presented static one-zone model suff'ers 
from limitations in both scenarios. In the lep tonic case, the 
model allows to estimate the average B -field (~ 25 jiG) and the 
total energy in accelerated electrons present in the shell of the 
SNR but fails to reproduce the observed X-ray and y-ray spec- 
tral slope. In the hadronic model, the high medium density re- 
quired to reproduce the observed TeV flux is hardly compatible 
with the hydrodynamics of the SNR. More detailed models using 
non-linear difl'use shock acceleration theory have been devel- 
oped (e.g. IZirakashvili & AharonianI 120 1 Ol: lEUison et al.l 120101) 
and would provide more accurate predictions than the simple 
model presented here. It should be noted that the considered 
model does not take into account evolution related to radiative 
cooling which could yield a steeper gamma-ray spectrum, in bet- 
ter agreement with the data. Also, the presented scenarios do 
not cover possible non-homogeneous surroundings such as wind 



Table 2. List of the parameters used for the spectral energy dis- 
tribution modeling presented in Fig. (TO] The spectral slope are 
fixed at 2.0 for the electron and the proton distribution. The den- 
sity of the ambient medium was set to 1 cm"^ in the case of the 
hadronic model. 



Model 








% 


B 




TeV 


TeV 


10^^ erg 


10^0 erg 


jiG 


leptonic 


18 




1.1 




25 


hadronic 


16 


100 


0.25 


2.0 


50 



bubble blown by the progenitor. Such detailed spectral and evo- 
lutionary modeling depends on many poorly known parameters 
and is therefore beyond the scope of the present discussion. 

Concerning the source HESS J1729-345, detected in the 
vicinity of the SNR, the presented multi- wavelength data do not 
provide a clear understanding of the nature of the object. The 
closest structures located near the y-ray emission are the HII 
region G353.381-0.114 (seen in the radio in Fig. |5l left) and a 
molecular gas clump observed in ^^CO (see Fig.[9l bottom) when 
integrating around a LSR velocity of -80 km/s (corresponding 
to near and far kinematic distances of ~6 and ~10 kpc respec- 
tively). If the y-ray source HESS J 1729 -345 is associated with 
those gas structures, it would therefore not be associated with 
the SNR HESS J173 1-347 thought to lie at a closer distance. 

6. Conclusion 

The newly discovered SNR HESS J173 1-347 exhibits a signif- 
icant shell morphology spatially resolved by H.E.S.S., similar 
to the one observed in radio. Together with RX J1713. 7-3946, 
RX J0852.0-4622 and SN 1006, HESS J173 1-347 is now the 
fourtljl TeV y-ray source to join this small but growing class. A 
lower limit to the distance of the SNR of 3.2 kpc was obtained 
by comparing the absorption derived from the X-rays and from 
HI and ^^CO observations. 

The multi-wavelength emission from the SNR, detected in 
radio. X-rays and y-rays, was combined in an SED to investigate 
the origin of the y-ray emission assuming that the broadband 
emission stems from the same region (one-zone model). While 
the measured fluxes can be accounted for in a purely leptonic 
model with a magnetic field of the order of 25 yuG, this simple 
model fails to reproduce the spectral shape of the X- and y-ray 
emission. A second model that assumes that the TeV emission is 
produced by hadronic processes is able to reproduce the spectral 
slopes in X- and y-rays at the cost of requiring that a large frac- 
tion of the kinetic energy of the explosion must be transferred to 
the accelerated protons and a high ambient medium density of 
^ ~1 cm"^ for J > 3.2 kpc. Moreover for such a density, the cor- 
responding shock speed of the SNR would be an order of mag- 
nitude lower than in other SNRs exhibiting bright synchrotron 
emission. 
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Fig. 10. Broadband SED for the sub-region of HESS J173 1 -347 
that is observed in X-rays (see Eig. [5] right, for region defini- 
tion). A purely leptonic (top) and a hadronic (bottom) scenario 
are shown and the corresponding parameters for both models are 
presented in Tabled The infrared (IR) seed photons energy den- 
sity and temperature were resp ectively set to 1 e V cm"^ and 40 
K following the model from Porte r et aP (|2008) for a galacto- 
centric radius of 4 kpc. The flux from the nearby Eermi source 
lEGL J1729.1 -3452c is represented and used as an upper limit. 
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